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Using Nonlinear Response to Estimate the Strength of an Elastic Network 
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Department of Physics, University of Houston, Houston, TX 77204 
The Institute of Fundamental Studies, Kandy 20000, Sri Lanka 

Disordered networks of fragile elastic elements have bee n proposed as a mod el of inner porous 
regions of large bones [Gunaratne et.al., cond-mat/0009221 , littp://xyz.lanl.gov|. It is shown that 
the ratio V of responses of such a network to static and periodic strain can be used to estimate 
its ultimate (or breaking) stress. Since bone fracture in older adults results from the weakening of 
porous bone, we discuss the possibility of using F as a non-invasive diagnostic of osteoporotic bone. 



Osteoporosis is a major socio-economic problem in an 
aging population [Q. Unfortunately, therapeutic agents 
which can prevent and even reverse osteoporosis often 
induce adverse side effects Thus, non-invasive di- 
agnostic tools to determine the necessity of therapeutic 
intervention are essential for effective management of os- 
teoporosis. Bone Mineral Density (BMD), or the effective 
bone density is the principal such investigative tool [^. 
Ultrasound transmission through bone [Q and geometri- 
cal characteristics of the inner porous region or trabecular 
architecture (TA) are being studied as complemen- 
tary diagnostics. 

In older adults, weakening of the TA is the princi- 
pal cause of increased propensity for bone fracture [Q. 
Analysis of models can complement mechanical studies 
of bone in aiding the identification of precursors of the 
weakening of a TA. In Ref. 1^, it was proposed that a sys- 
tem to model mechanical properties of a TA can be ob- 
tained by adapting a disordered network of fragile elastic 
elements . The model system includes potential energy 
contributions from elasticity and from changes in bond 
angles between adjacent springs. Furthermore, springs 
that are strained beyond (a predetermined value) e and 
bond angles that change more than 5 are assumed to 
fracture and are removed from the network. The iner- 
tia of the network is modeled by placing masses at the 
vertices. When in motion, each mass experiences a dis- 
sipative force proportional to its speed. Osteoporosis is 
modeled by random removal of a fraction v of springs 
from the network, and the bone density is estimated by 
the fraction of remaining links. Therapeutic regeneration 
is introduced by strengthening springs that experience 
large strain (as suggested by Wolff's law pO||). 

Numerical studies of the system show analogs of several 
mechanical properties of bone including, (1) an initially 
linear stress vs. strain curve that becomes nonlinear be- 
yond the fracture of elastic elements, (2) an exponen- 



tial reduction of the ultimate (or breaking) stress with 
decreasing BMD, and (3) a dramatic increase of bone 
strength following therapeutic regeneration . Together 
they support the conjecture that elastic networks are a 
suitable model of mechanical properties of bone. In this 
Letter we use results from a numerical study of the model 
to introduce a possible diagnostic tool for osteoporosis. 
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FIG. 1. The stress distributions on networks of size 60 x 60 
with (a) v — 10%, and (b) v = 30% representing "healthy" 
and "osteoporotic" bone respectively. For clarity only the 
compressed bonds are shown, and darker hues represent larger 
stresses. Notice that the "stress backbone" of (a) is dense, 
while that of (b) consists of a few coherent pathways. 

The primary difference between "healthy" and "weak" 
networks is the nature of stress propagation through 
them. Figure |^ provides a representation of stresses 
in two networks subjected to uniform compression 
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Figure 0(a) shows the effect on a "healthy" network 
(i/ = 10%) where large stresses supportmg propagation 
are seen to form a dense subset. In contrast, elastic el- 
ements supporting a "weak" network (i/ = 30%) form a 
few coherent pathways (Figure 0(b)). We refer to the set 
active in stress propagation as the "stress backbone" of 
a network JT2| . For a wide range of control parameters, 
it is seen to become finer with increasing v. 

It is easy to understand how the nature of the stress 
backbone is related to the stability of a bone. Loss of 
connectivity of a healthy TA (due to trauma) will have 
little effect on its stability, since many alternative path- 
ways are available for stress propagation. In contrast, 
fracture of a critical link (i.e., one belonging to the stress 
backbone) in a weak network will have a significant im- 
pact on the remaining stress pathways, likely inducing 
further fracture of elastic elements. Below, we discuss 
how these variations in stress backbones effect the non- 
linear response of networks to externally applied strain. 

Consider an elastic network from which a fraction v of 
elastic elements have been removed. It is first subjected 
to an adiabatically reached compression ■ Co is cho- 

sen well below the yield point so that there is no fracture 
of elastic elements. Under these conditions, the stress Tq 
needed to maintain the compression is proportional to 
Co |p^,p|, and the static susceptibility of the network is 
defined as xo = To/Co■ 
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FIG. 2. The power spectra of the dynamical susceptibility 
(with C,d = 0.001 and fio = 100) of several elastic networks 
normalized by the static susceptibility. The curves correspond 
to values of v of 0% (blue), 10% (green), 20% (red), 30% 
(aquamarine), and 40% (magenta). 

Next, this compressed network is subjected to an addi- 
tional periodic strain C(^) = Cp 6xp(ir2ot), and the force 



required for its implementation is denoted T,j{t). C^p is 
chosen to be small (Cp ^ Co)j and hence T^{t) can be 
assumed to be proportional to Cp- The dynamical suscep- 
tibility of the network is defined by Ti,{t) — Xv{t) ■ C{t)- 
The Fourier transform of Ty(i) is given by the convolu- 
tion T^{uj) — Xi,{u}) * C{(^)- Since C('^) — Cd'5(w -I- fio), it 
follows that xA(^) = Tu{uJ - ^o)/Cd 

Figure ^ shows the power spectra |Xi'('-^)l/xo for val- 
ues of 1^ of 0%, 10%, 20%, 30%, and 40%. Though both 
the static and dynamical susceptibilities reduce with ad- 
vancing "osteoporosis" (increasing v), Xi^i^) experiences 
a smaller reduction. This is possibly due to the forma- 
tion of additional (temporary) stress pathways when the 
network is subjected to periodic strain. 
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FIG. 3. The relationship between the ultimate (or break- 
ing) stress U{i') and T{i') for several elastic networks sub- 
jected to different compressions ^o. The symbols '-f ' (blue) 
and 'x' (green) represent distinct networks with identical con- 
trol parameters compressed by i^o = 4.0 and = 8.0 respec- 
tively. The symbols 'o' (red) and 'o' (aquamarine) represent 
two other network generated with different sets of control 
parameters, and compressed by = 5.0 and (^o = 3.0 re- 
spectively. For these networks the values of U{v) have been 
rescaled as described in the text. 

Variations of these power spectra can be quantified us- 
ing the ratio r(i/) defined by 

r-i(^) - - / Ix.MM^, (1) 

Xq J 

the range of the integral being a Nyquist frequency do- 
main |ll6|] . Figure ^ shows a remarkable relationship be- 
tween the ultimate stress U{i') and r(i^) for several net- 
works. The symbols and 'x' represent two distinct 
networks (constructed using different random seeds) with 
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identical control parameters which have been sub- 
jected to compressions of Co = 4.0 and Co — 8.0 units re- 
spectively. U (r) is seen lie on a common bi-linear curve. 

The behavior of U (F) for a third network is included 
to determine the effects of increasing the range of elas- 
tic moduli of the network The ultimate stress of a 
network is expected to reduce when the range of spring 
constants is increased (say, by a factor /). This is be- 
cause of the increase of the number of weak springs in 
the network. To compensate for this decrease, we rescale 
U{i') (heuristically) by/. Then, symbols 'o' representing 
the third network fall on the curve determined by the first 
pair of networks. U{T) will, of course, depend linearly on 
the mean elastic modulus, as has been confirmed. 

Finally, we include results from a fourth network to 
study variations of the fracture strain of elastic elements. 
A fixed-strain fracture criterion was used in the model; 
i.e., any spring that is strained by more than e, and and 
bond angle that is changed by more than 6 are removed 
from the network. This choice was motivated by obser- 
vations from recent mechanical studies which show that 
a trabecular architecture from a given location fails at a 
fixed level of strain, independent of the strength or elastic 
modulus of the bone [Q. However, bone samples from 
distinct locations exhibit different levels of fragility. For 
example samples from proximal rat tibia, human tibia, 
and human lumbar spine have been shown to fracture 
at strain levels of approximately 5%, 1% and 7% respec- 
tively 

The first three networks represented in Figure 3 in- 
cluded a common fracture criterion; specifically e = = 
5% and S = Sq = 0.1. Since F(^) is independent of e and 
S and U{i') can be expected increase with them, U{T) 
will depend on e and S. The symbols 'o' in Figure 3 rep- 
resent a fourth network |2^], where the ultimate stress 
U{i') has been rescaled by a factor (eo/ei), ei being the 
new value of the fracture strain. For parameters chosen, 
failure of elastic elements (as opposed to bond-breaking) 
was the dominant (though not exclusive) mode of frac- 
ture, justifying the use of this rescaling factor Once 
C/(j^) is rescaled, points representing all networks collapse 
to the same bi-linear curve p3[ |. 

The transition between the two linear segments of U (F) 
is accompanied by a qualitative change in the stress dis- 
tribution on compressed networks. Points on the right 
segment of each curve represent "healthy" networks; i.e., 
those with smaller values of ^. The histograms of stresses 
for these networks, an example of which is shown by the 
solid line of Figure 4, exhibit broad peaks. The presence 
of such peaks is consistent with the geometry of their 
stress backbones, see Figure 1(a). In contrast, points 



on the left segment of U{T) represent "weak" networks 
(i.e., larger v), whose stress histograms exhibit no broad 
peaks, and have exponential tails ||2^. This behavior of 
the histograms is consistent with the presence of a sparse 
stress backbones. The transition between the two linear 
segments of U{r) with increasing ly coincides with the 
elimination of the peak in the histogram. 
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FIG. 4. The solid line shows a histogram of stresses for a 
compressed network modeling a healthy TA = 10%). The 
dashed line shows the analog for a weak TA {v = 25%). Each 
histogram represents the average of 50 configurations with 
identical control parameters. The elimination of the broad 
peak of the histogram coincides with the transition between 
the two linear segments of U{T). Positive and negative values 
of the stress denote extended and compressed springs respec- 
tively. 

The ultimate stress [/ of a bone is the measure essen- 
tial for management of osteoporosis. Unfortunately, it 
is not accessible in- vivo (without breaking a bone!), and 
surrogates such as the bone density are used to estimate 
U. The BMD of a patient is compared with that of a 
sample population to determine if and when therapeu- 
tic interventions are necessary. However U is known to 
depend on other factors of bone such as the architec- 
ture of its TA and the "quality" of bone material. The 
significant variations introduced by these factors makes 
it difficult to identify individuals susceptible to fracture 
using measurements of BMD alone . 

These problems can be avoided if a characteristic that 
relates U to an intrinsic property (i.e., one that does not 
need to be compared to that of a population) of a bone 
is available. Since factors such as bone quality and archi- 
tecture of the porous medium effect both U and F it is 
conceivable that the relationship between them remains 
unchanged between patients. Numerical analyses of our 
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model system justify this expectation when U is rescaled 
by a location dependent factor (quantifying the stiffness 
and fracture strain of trabecular elements). 

Vibrational analysis has been used for in-vivo studies 
of bone strength and protocols are available to obtain 
measurements needed to evaluate F [^,^. Once the 
stiffness and fracture strain of different bone locations 
are tabulated, rescaling factors required for Figure 3 can 
be determined. Subsequently, the ratio F of responses of 
a bone to stationary and periodic strain can be used as 
a non-invasive diagnostic of bone strength. 
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ing out that nonlinear response is related to the stress 
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